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Abstract 

In this paper we study hypercomplex manifolds in four dimensions. Rather 
than using an approach based on differential forms, we develop a dual ap- 
proach using vector fields. The condition on these vector fields may then 
be interpreted as Lax equations, exhibiting the integrability properties of 
such manifolds. A number of different field equations for such hypercomplex 
manifolds are derived, one of which is in Cauchy-Kovaleskaya form which 
enables a formal general solution to be given. Various other properties of 
the field equations and their solutions are studied, such as their symmetry 
properties and the associated hierarchy of conservation laws. 
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1 Introduction 



The study of hyperKahler geometries has developed in two distinct directions. Starting 
with the work of Calabi |C| there has been a purely geometric vein, where such manifold are 
constructed and studied geometrically without reference to any defining set of field equa- 
tions. The second vein starts with such field equations - systems of differential equations, 
and uses solutions of such systems to construct the manifolds. 

To show how these two approaches are connected it is necessary to restrict one's at- 
tention to four dimensions, where the hyperKahler condition is equivalent to the existence 
of a metric with anti-self-dual Weyl tensor and vanishing Ricci tensor. Such metrics were 
shown by Penrose to have a corresponding twistor space, and conversely, that from such 
a twistor space one may reconstruct the metric. Such a structure appears at first sight to 
be special to four dimensions where there is the notion of self-duality, but it was shown 
in 1 5a | that hyperKahler metric have corresponding twistor space in any 4iV-dimensional 
space. 

With more recent work, initiated by Ward |W[ (see also ||MW|| ), these two threads 
may be seen to be intimately interwoven. The existence of a (suitable) twistor space 
indicates that one is dealing with an integrable system. Thus any set of field equations for a 
hyperKahler metric provides one with an example of a multidimensional integrable system. 
From a solution to this integrable system one may construct the associated twistor space 
(and vica- versa), whose properties may be studied, and hence properties of the metric, 
without recourse to the particular differential equation, whose precise form depends on the 
particular coordinate representation being used. 



The study of hypercomplex manifolds has had a similar history. In [FP|| Finley and 
Plebahski studied field equations which were derived from the following condition on the 
self-dual 2-form , 

dS* = qAS ! , i = 1,2,3, 



though without mention of the associated complex geometry. Similarly Boyer ||Bo|| studied 



the geometric aspects, but did not write down field equations for such metrics. In nei- 
ther case was the link with integrable systems made. The aim of this paper is to show 
that one may write down systems of differential equations, solution of which define hyper- 
complex manifolds (here we will restrict out attention to four dimensions, though many 
of the ideas will generalize to 4 TV-dimensions) . Since hypercomplex manifolds also have 
associated twistor spaces (in four dimensions the hypercomplex condition implies, though 
is not implied by, the anti-self-duality of the Weyl tensor | [Bo|| ) the field equations for these 



systems will be examples of multicomponent, multidimensional integrable systems. 

Since hyperKahler metrics are obviously Kahler they may be written in terms of a 
single functions, the Kahler potential Q : 

g = P dx i dx y , i,j = 12. 
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In four dimensions the hyperKahler conditions results in the differential equation 



^1^1 rj> 2 /j> 2 r£* 1 2 /ji 2 -f 



known as Plebahski's equation. In this form it is hard (though not impossible) to ap- 
ply ideas from the theory of integrable systems, which are best suited to evolutionary 
type equations. In ||G1|| the first author showed how, by performing a suitable Legendre 



transformation one may obtain an equation in evolutionary, or Cauchy-Kovaleskaya form, 
namely 

At = i>zt4>xy ~ IpyzTpxt (1-1) 



and the second author showed how to construct the associated integrable hierarchy, 
based on the study of the generalised symmetries of this equation. 

In this paper the following two-component generalization of (|1 . 1|) will be studied: 

9u = {gx,g} + {h,g z } , , 12 . 
h tt = {h x .g} + {h,h z } . 

From the purely integrable systems aspect, considered in section 3, one may view this 
system as resulting from dropping the volume preserving condition on the vector fields 
which appear in the Lax pair used to construct (|1 . 1|) . This approach will be followed 
in section 3. However, the system has a more geometric interpretation: solutions define 
hypercomplex metrics. This aspect will be considered in the next section. It should be 
pointed out that these hypercomplex structure are not the most general possible, but a 
particular subclass. The more general case will be considered in a future paper. 



2 Geometrical Description 

2.1 Hypercomplex geometry 

In this subsection, we wish to investigate a special subclass of four-dimensional hyper- 
complex structures. In particular, assume we are working on a four-manifold X, and that 
we have a local basis for the tangent space, in the form of a set of four linearly-independent 
vector fields, {ej : i — 1, . . . ,4}. We wish to study the geometrical structures which arise 
if we assume that these vector fields obey the Lie-bracket relations: 

[e 1 ,e 2 ] + [e 3 ,e 4 ] = 0, 

[ei,e 3 ] + [e 4) e 2 ] = 0, (2.3) 
[ei,e 4 ] + [e 2 ,e 3 ] = 0. 

We will show that in this case, the manifold X is (locally) hyper-complex. Recall that a 
manifold M of dimension n = 4m is hyper-complex if it admits three integrable complex 
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structures, I, J,K, which obey the quaternion multiplication relations: 

I 2 = J 2 = K 2 = -ld TxM , 
IoJ = K, JoK = I, KoI = J. 

Such structures imply that the bundle of linear frames L(M) reduces from a GL(4m,M.) 
bundle to a GL(m, H) bundle. In the special case of four dimensions, the frame bundle 
reduces to a GL(1,H) = W = R + x SU(2) bundle. Since 577(2) is a subgroup of SO(A), 
this means that in four dimensions, a hypercomplex structure automatically defines a 
conformal structure. In particular, there is a metric, unique up to a scale, with respect to 
which all the complex structures are Hermitian: 

g(IX,IY) = g(JX, JY) = g(KX,KY) = g(X, Y), 

where X, Y are arbitrary sections of TM. The structures I, J, K along with a representative 
metric in the conformal structure, g , define a hyper-hermitian structure. Each metric 
in this conformal structure has anti-self-dual Weyl tensor, with respect to the canonical 
orientation defined by any of the complex structures. 

In the case we wish to study, we have a set of vector fields which satisfy the relations 
of equation Q2.3Q . If we define the dual basis {e 1 } for T*X, then we will show that the 
vector fields define three integrable complex structures which obey the relations given in 
equation fl2.4| ) and that the metrics with respect to which all these structures are hermitian 
are conformal to the metric: 

g = e 1 <g> e 1 + e 2 <g> e 2 + e 3 <g> e 3 + e 4 <g> e 4 . 

The plan of the section is as follows. In the next section, we prove the assertions made 
above concerning the existence of hyper-complex structures, and a compatible conformal 
structure. Next the relation between this approach, based on vector fields, and the more 
usual approach based on forms is given. The vector field approach is useful for two reasons: 
firstly it enables field equations to be derived easily, and secondly it makes the connection 
with integrable systems more transparent. We then consider various coordinate versions 
of these equations. 

Theorem 2.1 On a four- dimensional manifold X , if there exist four linearly-independent, 
non-vanishing vector fields {ej : i = 1,2,3,4} on a manifold X which obey the Lie Bracket 
relations 

[ei,e 2 ] + [e 3 ,e 4 ] = 0, 

[ ei ,e 3 ] + [e 4 ,e 2 ] = 0, (2.4) 
[ei,e 4 ] + [e 2 ,e 3 ] = 0, 

then the manifold is locally hypercomplex. Moreover, there exists a unique conformal struc- 
ture, defined by the representative metric: 

g = e 1 <g) e 1 + e 2 <g> e 2 + e 3 <g> e 3 + e 4 <g> e 4 . (2.5) 

with respect to which all of the complex structures are Hermitian. 
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Proof Consider maps (I, J, K) on TX defined by 

I(ei) = +e 2 , I(e 2 ) = -e 1; I(e 3 ) = +e 4 , I(e 4 ) = -e 3 , 
J(ei) = +e 3 , J(e 2 ) = -e 4 , J(e 3 ) = -e x , J(e 4 ) = +e 2 , 
K(ei) = +e 4 , K(e 2 ) = +e 3 , K(e 3 ) = -e 2 , K(e 4 ) = -e^ 

To impose that these structures are integrable, we consider the complexification of the 
tangent space T C X = TX <g) C. Consider first the structure I. The almost complex 
structure leads to a direct sum decomposition of T C X as T^ 1 ' ) © T^ ' 1 ) where v e T^ 1 ' ) if 
I(v) = iv, and v 6 if l(y) = —iv. The complex structure is integrable iff the space 

T^- 1 ' ' is closed under the Lie bracket. In the case of the structure I, T*- 1,0 - 1 is spanned by 
the vector fields {e 4 + ie 2 , e 3 + ze 4 }, and the only non-trivial Lie bracket we must consider 
is 

[ex + ie 2 , e 3 + ze 4 ] = [e 4 , e 3 ] + [e 4 , e 2 ] + i ([e x , e 4 ] + [e 2 , e 3 ]) . 

For the right hand side to lie in T^ 1 ' ' implies a number of algebraic relations which, together 
with the analogous equations for J and K and the quaternionic relations, implies that 

[ei, e 2 ] + [e 3 , e 4 ] = -A 2 e 1 + A x e 2 - A A e 3 + A 3 e 4 , 
[ei, e 3 ] + [e 4 , e 2 ] = -A^ + A A e 2 + Aie 3 - A 2 e 4 , 
[ex, e 4 ] + [e 2 , e 3 ] = -A A e x - A 3 e 2 + A 2 e 3 + A x e A 

for some set of functions {Ai, A 2 , A 3 , A 4 } . Thus by construction the corresponding metric 
( p.5|) is hypercomplex. It will be useful in what follows to combine these functions into a 
one form A = A^ and regard this as a connection on X . 

Unlike the conditions for a metric to be hyperKahler, the hypercomplex conditions are 
invariant under conformal changes of the metric. If g — > e 2A g then the corresponding 
transformation for the connection A is 

A -> A + 2dA . 

The subclass of hypercomplex structures which we are examining in this paper are defined 
by the conformal invariant condition dA = 0, i.e. the connection defined by A is flat. 
Thus locally one may define, for this subclass of hypercomplex structures, the conformal 
factor so that A = . This fixes the conformal structure. 

It is a straightforward exercise using the explicit form of the complex structures given 
above to show that metric g given above is Hermitian with respect to each of the complex 
structures, and that it is the unique symmetric tensor (up to a rescaling) with this property. 

□ 

The hypercomplex condition determines a metric up to a conformal factor. The zero 
curvature condition on the connection A may be used to fix the conformal factor. This 
determines the metric uniquely, up to trivial transformations. In terms of a null tetrad in 
which 

g = e 1 0s e 2 + e 3 ® s e 4 
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the conditions ( |2.4| ) become 



[ei,e 2 ] + [e 3 ,e 4 ] = 0, 
[ei,e 3 ] = 0, 
[e 2 ,e 4 ] = 0. 

This form will be used in later sections and is also used in the following example. 
Example Consider the following vector fields 



(2.6) 



e-i 
e 4 



U 



e 2 
e 3 



[1 + \w\ 2 )d w + zwdz, 
1 + \z\ 2 )dj + zwdm . 



It is easy to verify that these satisfy the conditions ( |2.6| ). Hence, by the above theorem, 
the corresponding Hermitian hypercomplex metric on C 2 is conformal to 

g = (1 + \z\ 2 )dw dw + (1 + |z| 2 )<iz dz — zwdz — zwdz dw . 

Note, this metric is locally conformal to the Fubini-Study metric on CP 2 . However this 
construction does not extend from C 2 to CP 2 . If it did it would contradict Boyer's [|Bo|| 
classification of compact hyperHermitian manifolds. Other examples, with tri-holomorphic 
Killing vectors, have been constructed in ||CTV| , |GT . 



For the metric to be hyperKahler the corresponding Kahler forms defined by 



n T (x,Y) 
nj(x,Y) 
n K (x,Y) 



g(IX,Y), 
S{JX,Y), 
S(KX,Y), 



must be closed, or equivalently, that V/ = VJ = VK = . It was shown in [MJN | that 
these conditions are equivalent to the vector fields being volume preserving, that is 

C-.uj = 



where u) is some volume form. 



2.2 Dual Description 

The starting point for the study of four dimensional hypercomplex manifolds has tradi- 
tionally been the equation 

dE^aAS', i = 1,2,3 (2.7) 

where the IT are self-dual two forms on the manifold X . In this paper we have so far used 
a dual description, using vector fields rather than forms. This subsection is intended to 
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bridge the gap between these two approaches. It is first necessary to fix some notation. 
The connection one-forms are defined by 

de* + rv A e j = , 

so, in components, Y l j = r* -^e* . The antisymmetric parts of P - fe are related to the 
structure functions defined by the Lie bracket [ey,efe] = c- fc *ej by 

This formula enable one to connect these two approaches. 



Proposition 2.1 The connection between equations t \2.3j ) and ( \2. 7y is given by the for- 
mulae 



ol = A - x 



where 



A = Aie 1 , 



X = c w V. 



Proof Consider the self-dual two-form S 1 = e 1 A e 2 + e 3 A e 4 . Then 



dS 1 = d(e l A e 2 + e 3 A e 



-r^e" A e b A e 2 + T 2 [ab] e a A e b A e 1 — T% b] e a A e b A e 4 + T 4 [ab] e a A e b A e 3 
-ic o6 V A e b A e 2 + ic o6 V A e 6 A e 1 — I c<lk V A e b A e 4 + ± Cfl6 V A e 6 A e 3 

+6 1 A e 2 A (c afc a e fe + A 3 e 3 + A^) + e 3 A e 4 A (c ab a e b + A.e 1 + A 2 e 2 ) 
(e 1 A e 2 + e 3 A e 4 ) A (A - X ) , 
(A - X ) A S 1 . 



The manipulations for the remaining self-dual forms are identical 

□ 

The hypercomplex condition on the metric is invariant under conformal changes of the 
metric. Thus one needs to calculate the transformation properties of these forms under 
such a change. 
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Lemma 2.1 Under the conformal change g — > e g the forms a , A , x transform as 

a — > a + 2dA , 
A -> A-dA, 
X -> X-3dA. 

Proof The proof of this is entirely straightforward, following directly from the definition of 
the various forms. It is interesting to note that from these one may construct two linearly 
independent conformally invariant two forms. 

□ 

Locally, if da = then a conformal factor may be found so the resulting metric is hy- 
perKahler. The metrics studied in this paper come from the conformally invariant condition 
dA = . This leaves another conformally invariant condition d\ = , the significance of 
which will be investigated elsewhere. 

2.3 Local coordinate representations 

The aim of this section is to give some local coordinate descriptions of the equations (|2.4j ). 
The main point that we make use of is that the equations ( |2.4j ) may be interpreted as the 
integrability conditions for null planes in the complexified tangent space of the manifold, 
and that we can introduce coordinates which naturally describe these surfaces. 

With this approach in mind, we consider the complexified tangent bundle, T C X = 
TX eg) C. From the vector fields {ej}, we define the complex basis for T C X: 



u 


= ei + ie 2 


v 


= e l - ie 2 


w 


= e 3 + ze 4 


x 


= e 3 - ze 4 



In the case where we are considering the complexification of a real Riemannian structure, 
we would impose the reality condition that u and v would be related by complex conju- 
gation, as would w and x. More generally, we require the existence of a real structure on 
T C X to define real slices of different signatures. Since we will be interested in metrics of 
both Riemannian and ultra-hyperbolic signatures, we will generally assume, for the mo- 
ment, that the vector fields u, v, w,x are complex, and simply impose the relevant reality 
conditions later. 

In terms of the vector fields u, v, w, x, the relations ( |2.4| ) take the form: 

[u, w] = 0, 
[u,v] + [w,x] = 0, 
[v,x] =0. 
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The first of these equations states that at each point (on a suitable local neighbouhood of) 
X, the plane spanned by the vector fields u and w is integrable. It follows that we can 
introduce coordinates (t, y) on each of these surfaces with 

u = d t , w = d y . 

If we now consider the space with t and y held constant, the last of equations ( |3.10 ) 
tells us that this space forms an integrable plane as well. As such, we can introduce 
coordinates (x, z) on these planes, and the coordinates (t, x, y, z) should give a suitable 
coordinate system on some local region of X. Note that the coordinates ( not 
fully determined by these conditions. The different choices of coordinates correspond to 
different coordinate expressions for the vector fields v and x, and in fact there are several 
geometrically distinct coordinate systems that we wish to investigate: 



2.3.1 Case I 

The first case we consider is where we take the coordinates (x, z) to be null. This is the 
analogue of the usual complex coordinate description of Kahler metrics. Using the second 
of equations (|3.1U| ), we introduce functions a and b and let v and x take the form 

v = a y d x -b y d z , x = -a t d x + b t d z . 

This is the generalisation of the form used in [|CMN| for the case of half-flat metrics. The 
functions a and b must now satisfy the equation 

{a, a x } = {b, a z }, {a, b x } = {b, b z }, (2.8) 

where we have defined the Poisson Bracket by 

{f,g} = ft9 y - f y 9t- 



Equations (|2.8|) are a generalisation of the First Heavenly equation which describes half- 
flat metrics JFJ. To see this, we note that the vectors in ( |2.4|) define a metric which is 
conformal to a half flat metric if the vectors are divergence free with respect to some 
volume element uj | MN |. Taking uj = dt A dx A dy A dz, we find that there exists a function 



Q such that a = Q z ,b = Q y . We can therefore integrate equations fl2.8|) once and rescale 
our coordinates (y, z) such that Q satisfies the equation 

{n y ,n z } = 1, 

which is the First Heavenly equation |FT|. 

In this case, we may reconstruct the metric, g, which takes the form: 

g = 4 (b t a y — dtby)' 1 [dt ® (a t dz + b t dx) + dy ® (a y dz + b y dx)\ . 

We see that in this form, the coordinates (t,x,y,z) are all null, and are tailored to the 
geometry of the spheres worth of integrable null planes through each point of the space. 

The metric is Riemannian if we assume that the functions a and b are real, and that 
the coordinates obey some reality condition. If the functions and the coordinates are real, 
the metric will be of ultra-hyperbolic signature. 
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2.3.2 Case II 



Here we introduce functions (p and ip and take 

v = d x + 4> y d t - tpydy, x = d z - <p t dt + ifitdy. 
The functions <fi and if) must now satisfy the equations 

4>tx + 4>yz + {<fit, 4>} + {i>, 4>y} = o, 

ifitx + Ipyz + {ik, 4>] + {ifi, Ipy} = 0. 

In terms of these coordinates and functions, the metric is 

g = 4 [(dt + (fi t dz — (fiydx) ® dx + (dy — ifj t dz + ip y dx) <g> dz] . 

In this case, the coordinates (t, y) are again null labelling some of the null planes in the 
space. The coordinates (x, z) are not null, however, and label the "rate of change" of the 
null planes (see [|NP'1| for more on the geometrical interpretation of these coordinates.) 

These equations have been previously investigated in connection with anti-self-dual 
structures |[FP|| , and are a direct generalisation of the Second Heavenly Equation for half- 



flat metrics IF! 



2.3.3 Case III 

Finally, we consider the analogue of the expansion used in to reduce the half-flat case to 
evolution form [pl|| . We therefore introduce functions g and h such that 

v = d t + g y d x - h y d z , x = -g t d x + h t d z . 

We find that the functions / and g must satisfy the equations 

gtt + {gx,g} + {h,g z } = 0, h tt + {h x ,g} + {h, h e } = 0. (2.9) 

In this coordinate system, the metric takes the local form 

g = 4A _1 [dt (g> (gtdz + h t dx) + dy <g> (g y dz + h y dx) — A" 1 (gtdz + h t dx) 2 ] 

with A = (h t g y — gth y ) . In this system, the coordinates (t, y) still span null planes, how- 
ever the geometrical interpretation of the (x, z) coordinates is not particularly clear. The 
main advantage of this coordinate system, however, is that the equations of motion are in 
evolution form. This form of the equations is therefore the natural starting point for the 
study of symmetry algebras ||G1|| , and the associated integrable hierarchy We shall 



therefore concentrate on this form of the equations from now on. 
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3 Integrable Description 



In this section the system ( |1.2|) will be studied, viewing it as an integrable system and hence 
applying various known results from the theory of integrable systems to it. In particular 
a Lax pair will be given for the system, a hierarchy of conservation laws constructed and 
the Lie-point symmetry structure calculated. 



3.1 Lax Pair 

A characteristic feature of an integrable system is the ability to express it as the com- 
patibility condition for an otherwise overdetermined linear system. To obtain the system 



1.2 ) in such a way consider the following vector fields on complexified tangent bundle, 



T C X = TX ® C of a four manifold X 



u = d t , w = d y , 

v = d t + g y d x - h y d z , x = -g t d x + h t d z . 



With these define the new vector fields 



C = u — Ax , 
C\ = w + Av 



where A G CP 1 is an auxiliary parameter. The compatibility conditions for the otherwise 
overdetermined linear system — = results in the following system of equations: 



[u, w] = 0, 
[u, v] + [w,x] = 0, 
[v,x] =0. 

With the explicit vector fields given above these reduce to the system ( |1.2j ): 

9tt = {g x , g} + {h, g z } , ,^iq\ 
h a = {h x ,g} + {h, h z } , 

where, for convenience, the Poisson bracket {/, g} = f y g z — f z g y has been used. In the 
previous section the geometry underlying this construction was given; solutions define a 
hypercomplex metric on the manifold X . Here we just consider the system as an example 
of a four dimensional integrable system and study it thus. 

One interesting reduction of this system is to impose the condition g y = h z . One may 
solve this constraint by introducing a function if; such that h = ip y , g = ip z . With this 
it is possible to integrate ( |L2| ) and obtain a single evolution equation fll.ip . With this 
constraint the basic vector fields ei = (u , v , w , x) become volume preserving, that is 

C ei u> = 
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where uj is the volume form uj = uAvAwAx. In this case the metric is hyperKahler 
rather than hypercomplex, with ip being related to the Kahler potential by a Legendre 
transformation. Conversely, one may regard the system (|1.2|) as a generalization of 
when one relaxes the volume preserving condition on the vector fields in the associated 
Lax pair. 



3.2 Formal Solutions 



The system ( |1 . 1|) is in Cauchy-Kovaleskaya form, so their formal solution may be written 



as a power series in the t-variable: 



h = ^2h n (x,y,z)t n 

71=0 
OO 

9 = ^2g n (x,y,z)t n 



n=0 

and the differential equations reduce to recursion relations between the coefficients h n , g n . 
Thus a formal solution may be derived from the coefficients (go ,ho,gi,hx),ox equivalently, 
in terms of the initial data (g ,h,g t , h t )\ t=0 on the t — hypersurface in X . This shows 
that the general solution depends on four arbitrary functions of three variables. 

One simple, but explicit, solution may be obtained from the ansatz 



g = ty + G(t,x,z) , 
h = t. 

With this the nonlinearites in (|l.lj) disappear and one is left with a linear equation for G , 
which after a simple change of variable, is just the three dimensional Laplace equation. 
Other simple solutions may be obtained by taking known hypercomplex metrics and reex- 
pressing them in the above form. Some examples of solutions obtained in this way will be 
given later. 

3.3 Symmetry Structure 

While symmetry techniques may be applied to any system of differential equations, the 
Lie-point symmetries of integrable systems have a particularly rich structure compared 
to non-integrable systems. Indeed, possible integrable systems may often be indentified 
by an increase in the dimension of the Lie-algebra of symmetries, as compared to nearby 
non-integrable systems. 
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Let ) and u = (-u 1 , . . . , u q ) be sets of independent and dependent 

variables, and consider a set of differential equations of degree k, given by 

A*(x,u (fc) ) = 0, i = l,...,m. 

Lie-point symmetries are generated by the vector field 

p d q d 

where the coefficients are determined by the criterion 

pr«v(A)| A=o = 

where pr^v is the fc-th prolongation of the vector field v . These ideas, and notation, are 
standard, see fOR . 



To apply such a procedure to the system ( |1 . 2| ) it is convenient to convert the system 
from a second order system in two independent variables to a first order system in four 
independent variables by introducing potentials and integrability conditions. Therefore let: 

A = g t ,B = g y ,C = h t , D = h y , 

in which ( Q[) becomes 

A y -B t = 0, 

Cy - D t = 0, , , 

A t + A X B-B X A + CB Z -DA Z = 0, > 
C t + C X B - D X A + CD Z -C Z D = 0. 

Since these are first order the calculation of the first prolongation of the vector field is easy, 
the evolutionary form of the system also giving a distinguished variable t to eliminate in 
the course of the calculations. 

The result is that we get five families of symmetries, the vector fields which generate 
these families being (where A; is a constant, cj), if) are functions of the coordinate y, and 
ai, 0,2 are functions of coordinates (x, z))\ 



v 1 ^] = k(-td t + 2Ad A + Bd B + 2Cd c + Dd D ) 

yd B + C(f)yC 



v 2 [0] = - ( f ) d t + A(f> y d B + C(f> y d D , 



V 3 [^] = -tlP y d t -4>dy+?Py(Ad A + Bd B +Cdc+Dd D )+t?Py y (Ad B + Cd D ), 

v 4 [ai] = -axd x - {a Xx A - a lz C) d A - {ax x B - a lz D) d B , 
v 5 [a 2 ] = -a 2 d z - (a 2z C -a 2x A)d c - (a 2z D -a 2x B)d D . 

The first of these generators simply generates a scaling symmetry of the equations. There- 
fore, if A, B, C, D constitute a solution of the equations, then so do A, B, C, D defined 
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by 

A(t,x,y,z) = e 2k A(e k t,x,y, z), 

B(t,x,y,z) = e k B(e k t,x,y,z), 

C(t,x,y,z) = e 2k C(e k t,x,y,z), 

D(t,x,y,z) = e k D(e k t,x,y, z). 

Similarly, the second generator generates a translation in the t coordinate along with a 
redefinition of fields, so that in this case 

A(t,x,y,z) = A(t + <f>,x,y,z), 

B(t,x,y,z) = B(t + <f>,x,y,z) + </> y A(t + (/>,x,y,z), 

C(t,x,y,z) = C(t + (j>,x,y,z), 

D(t,x,y,z) = D(t + <j),x,y,z) + <p y C(t + <j),x,y,z). 

These are the only symmetries that it is possible to generally exponentiate explicitly. The 
non-zero commutators for this Lie algebra are: 

[v^v 2 ^]] = V 2 [fc0], 
[v 3 [^],V 3 [x]] = V 3 bcip y - l/lXy], 

[v 4 [a],v 4 [6]j = v 4 [a x b-ab x ], 
[v 4 [a],v 5 [6]] = v 4 [a z b]-v 5 [ab x ], 
[v 5 [a],v 5 [6]j = v 5 [a z b-ab z ], 

where k, I are arbitrary constants, -0, are arbitrary functions of y, and a, b are arbitrary 
functions of x and z. 

From the structure of these commutators one may decompose the Lie algebra L gener- 
ated by these vector fields into a direct sum L = L 1 © L 2 , where 

L 2 = {v 4 [a],v 5 [6]}. 

The sub-Lie-algebra L\ decomposes further as a semidirect product L\ = S > R where 

S = {v 3 ^}, 

R = {v^/c], v 2 [0]} . 

The subalgebra L 2 is isomorphic to vector fields on a 2-dimensional surface and corre- 
spond to coordinate transformation in the x , z-variables. This is to be expected, since in 
the hyper-Kahler case the vectors are all divergence free and the symmetries turn out to 
be related to symplectic diffeomorphisms of two dimensional planes. In our case, we have 
simply dropped the divergence-free condition from the vector fields, and the symmetry 
group becomes related to the larger group of diffeomorphisms, since there is no natural 
symplectic structure any more. Similarly the vector fields v 2 generate coordinate trans- 
formations and the vector fields v 1 generate constant rescalings of the metric. The only 
vector fields which generate genuinely new metrics are those in S . 
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3.4 Conservation Laws 



In this section a hierarchy of conservation laws of the form 

<rv^t } = o, n = o,i,. 



will be constructed. This expression is clearly covariant, but in the calculations it will 
be necessary to use a particular form of the metric and the associated field equations. 
Explicitly we consider metrics of the form 



g = 2A 2 (g) [a t dz + b t dx) + dy <g> [a y dz + b y dx)] 



(3.12) 



where A = (a t b y — a y b t ) and with a and b being solutions of the field equations (|2.8|) . The 
conformal fact in ( |3.12j ) has been fixed so that det = 1 ; such a fixing does not change 
the hypercomplex or Hermitian properties of the metric. One obvious extension of these 
results would be to introduce the notion of a conformally invariant conservation law. 

The starting point of this construction, a generalization of a procedure first applied to 
nonlinear cx-models |[B1ZZ| , CP], is the solution \& to the Lax pair 



[\d t + a t d x -b t d z ]V = 
[Xd y + a y d x - b x d z \ * = 



A e CP 1 . 



(3.13) 



Expanding $ as a power series in A , i.e. as \1/ = Yl^Lo ^ n ^n an d equating coefficients 
yields the following equations for the \I/o _ term: 



[a t d x -b t d z }^ = 0, 
[a y d x - b y d z ) # = 0, 



and the recursion relations 



dt^n 



-a t d x + b t d z )^ n+ i , 
-a y d x + b y d z )^f n+1 . 



(3.14) 



(3.15) 



The first set of equations imply that \l/ = ^o(^l/) an d so we take the seed solution to be 

V y 

(here we have assembled two independent solutions into a vector). This seed solution will 
generate, via the recursion relations ( |3.15 ) the full solution to the Lax pair ( 3.13| ). This 
function defines the twistor surfaces in the corresponding twistor space. Another family 
of conservation laws may be obtained starting from the expansion \1/ = Yl^=o ^~ n ^n and 
the relationship between \1/ and \I/ on the equator of CP 1 defines the twistor space, via a 



patching construction ||NPT| , |S2 
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Proposition 3.1 The currents jp, defined by 

A n) = 0, jf } = A^ n+1 , 

are conserved. 

Proof With the particular metric ( |3.12| ) 

g^nJ^ = d t [-ayd z y n+1 + b y d z y n+1 ]+dy[+atd z y n+1 -b t d z y n+1 ], 
= d t [dS> n ]-d y [d t V n }, 
= 0. 

This proof uses the condition det = 1 , so the Christoffel symbols T% = . 

□ 



4 Comments 

Underlying the integrability of the multidimensional systems presented here is the existence 
of a twistor space. This paper has, though, only concentrated on the field equations and 
the associated Lax pairs with little mention of the properties of the corresponding twistor 
space - in terms of a double fibration 

{Lax pair} 

/ \ 
{hypercomplex manifold} {twistor space} 

we have said little about about the structure of the right hand side. Such twistor space have 
the special property that they fibre over CP 1 ra, unlike those for more general ant i- self- 



dual Weyl spaces or scalar flat Kahler spaces. This is manifested in the simple A-dependence 
in the Lax pairs for hypercomplex manifolds - the Lax pairs for scalar flat Kahler i pr'aj 



and general anti-self-dual Weyl spaces [|G2| , |MW|| involve c^-terms. The hypercomplex 



manifolds studied here come from the conformally invariant condition dA = , and so one 
would expect the corresponding twistor space to exhibit certain extra properties. General 
hypercomplex manifolds (without this condition) will be studied in the sequel to this paper, 
this also containing the connection between the approach developed here and the Obata 
connection. 

One characteristic feature of integrable systems is the existence of an associated hierar- 
chy. Such hierarchies may be constructed by studying the generalized symmetry structure 
of the original systems of equations [|OJ. Such hierarchies have been constructed for hy- 
perKahler metrics in |52|. It remains to see how such ideas may be extended to the 



hypercomplex systems studied here. 
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Appendix 

While many of the ideas in this paper will generalize to 4iV-dimensional manifolds, when 
N = 1 another generalization is possible. Underlying the integrability of the structures 
studied in this paper is the existence of a suitable twistor space, and hypercomplex mani- 
folds automatically have such twistor spaces. However in four dimensions the existence of 
a twistor space follows from the Weyl tensor being anti-self-dual, and the hypercomplex 
condition here implies, not is implied by, this condition. Thus a possible generalization 
is to study metrics with anti-self-dual Weyl tensor and which are also scalar flat. Such a 
system has an other reduction to scalar flat Kahler metrics. These different systems and 
their interconnections are summarized in the following diagram: 



a.s.d. Weyl 
(R = 0) 

hypercomplex 1 J Kahler 

(R = 0) J \ R=0 

\ / 
HyperKahler 

Ricci flat 

The conditions (R = 0) in brackets indicate how the conformal factor for otherwise con- 
formally invariant conditions have been fixed. 

An analogous system of equations to ( |1.1|) for the scalar-flat, anti-self-dual- Weyl systems 
is given by 

9tt + (e-^{g,h}) z = e-^{g,h z -g y }, 
h tt +(e-^{g,h}) y = e-^{h,h z -g y }, 
{e + % t + {g^ y }-{h^ z } = 

(these certainly imply the geometric conditions though whether they are implied by them 
is unclear). The corresponding metric is given by 



+2dy< h t dt + h x dx —[h t dy + g t dz] > + 2dz< g t dt + g x dx -^-[h t dy + g t dz] } . 

where A = h t g x — gth x . The two reductions above are easy to see from this system: 
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when ip = , this system reduces to 



9u = {h,g z } - {g,g y } , 
h tt = {h,h z } - {g,h y } , 

that is, to the hypercomplex systems studied in the main body of this paper. The 
further reduction h = 9 y , g = 9 Z reduces this down to the hyperKahler equation 

(0)- 

• imposing the Kahler condition on this system gives h = 9 y , g = 9 Z and the first set 
of equations simplify to 



9 u + e-^{9 z ,9 y } = 0, 
(e + %t + {9 z ,ipy}-{9 y ,ip z } = 0. 

These are the analogues of the well-known scalar flat Kahler equations , written 
in evolutionary form. It is easy to find solutions, such as the one which gives the 
Burns metric. The further reduction ip = reduces this down to the hyperKahler 
equation ( |1.1|) . 

One interesting class of solutions to all these systems comes from imposing an SU(2) 
symmetry on the metrics. Such metrics are often referred to as Bianch IX metrics. This 
symmetry reduces the field equations from partial differential equations down to systems 
of coupled ordinary differential equations which may be integrated directly. These ideas 
may also be applied to other Bianchi metrics. 
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